Spectral function of a Luttinger liquid coupled to phonons and angle-resolved 
photoemission measurements in the cuprate superconductors 
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We compute the finite-temperature single-particle spectral function of a one-dimensional Luttinger 
liquid coupled to an optical phonon band. The calculation is performed exactly for the case in which 
electron-phonon coupling is purely forward scattering. We extend the results to include backward 
scattering with a renormalization group treatment. The dispersion contains a change in velocity at 
the phonon energy, qualitatively similar to the case of electron-phonon coupling in a Fermi liquid. 
If the backward scattering part of the the electron-phonon interaction is not too strong compared to 
the forward scattering part, coupling to phonons also produces a pronounced peak in the spectral 
function at low energies. The calculated spectral function is remarkably similar to the angle- resolved 
photoemission spectra of the high-temperature superconductors, including the apparent presence of 
"nodal quasiparticles," the presence of a "kink" in the dispersion, and the non-Fermi-liquid frequency 
and temperature dependencies. Although a microscopic justification has not been established for 
treating the electronic dynamics of the cuprates as quasi-one-dimensional, at the very least we take 
the quality of the comparison as evidence of the non-Fermi-liquid character of the measured spectra. 



I. INTRODUCTION 

There is increasing, although controversial, experimen- 
tal evidence that, in cuprate high-temperature supercon- 
ductors, a strong coupling between electrons and optical 
phonons produces observable features in the single-hole 
spectral function A(k, E) measured by angle-resolved 
photoemission spectroscopy (ARPES) iSiii In attempting 
to understand the implications of these measurements, 
the experiments have been previously analyzed in the 
context of the theory of the electron-phonon (el-ph) cou- 
pling in a Fermi liquid. However, the observed spectral 
function exhibits a manifestly non- Fermi-liquid frequency 
and temperature dependence so the justification for 
this mode of analysis is not clear. 

In the present paper, we compute and analyze the 
finite-temperature spectral function of a spinful, one- 
dimensional (ID) Luttinger liquid (LL) coupled to a 
dispersionless optical phonon band (Einstein phonon). 
Such a spectral function was computed previously in 
Ref. H but only at T — and in the absence of both 
electron-electron (el-el) interactions and el-ph backscat- 
tering. The LL is a quantum critical point, so the spec- 
tral function is a scaling function, and should be com- 
puted at finite temperature; the zero-temperature spec- 
tral function reveals only the high-frequency behavior of 
the scaling function. It is also essential to include el-el 
interactions, as in real materials they are typically much 
stronger than the el-ph interactions. Moreover, strong el- 
el interactions make qualitative changes to the "appear- 
ance" of the spectral function. The effects of phonon- 
assisted backward scattering are also important to con- 
sider, as we have done below. 

Our motivation for this study is threefold. (1) There 
are a host of interesting quasi- ID materials j2*i2i some of 
which are amenable to ARPES studiesjii*i2ii£ to which 
this analysis may be directly applicable. (2) The LL is 
the theoretically best understood example of a non Fermi 



liquid-just as many aspects of the Fermi liquid state are 
robust, independent of details of material and even di- 
mensionality, we may hope that some features of non- 
Fermi-liquids are similarly generic, at least within classes 
of non-Fermi-liquids. In this case, the LL may serve as a 
paradigmatic model for a broader class of systems. (3) It 
may be the case that the quasi-2D cuprates have signif- 
icant, self-organized "stripe" structures^ which render 
them locally quasi-lD, in which case it may be possi- 
ble to directly compare the results obtained here with 
experimental in the cuprates and other highly correlated 
materials. 15 

Let us start with a qualitative summary of the solu- 
tion of the el-ph problem in a Fermi liquid^ at zero tem- 
perature. Consider the case of an optical phonon with 
frequency ujq and dimensionless el-ph coupling A', which 
is not too large. For \E\ > u>q, the effect of the el-ph 
coupling on A(k, E) is an ^-independent broadening of 
the quasiparticle peak. For \E\ <C loq, the phonons can 
be integrated out to produce new effective interactions in 
the Fermi liquid: the largest effect is a renormalization of 
the Fermi velocity, Vf — ► v F = + A') <«f, while 

the most dramatic effect is the weak effective attraction 
produced between low-energy quasiparticles, which can 
lead to a superconducting instability of the Fermi liquid 
state. 

We have computed the single-hole spectral function for 
a LL coupled to optical phonons for the case when the 
el-ph coupling does not produce a spin gap. We find that 
the same words describe the effects of the el-ph coupling 
as in the Fermi liquid case, but the results look quite 
different because the unperturbed LL spectral function 
is not a simple Lorentzian. Specifically, if we let {g} 
represent the set of coupling constants which define the 
LL (i.e. the charge and spin velocities v c and v s , and the 
corresponding Luttinger exponents K c and K s ), then our 
result can be summarized by 
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A(k,E) 



A LL {k,E-{g}) + 
A LL (k,E;{g*}) 



for \E\ > ujq, 
for \E\ < w , 



(1) 



and there is a smooth crossover between the two limits 
when \E\ ~ u>n- Here, is the spectral function of the 
pure LL (i.e. in the absence of el-ph coupling), the el- 
lipsis represents perturbative corrections to the spectral 
function which can be ignored as long as the cl-ph cou- 
pling is not too big or if \E\ is high enough, and {g*} are 
renormalized coupling constants obtained by integrating 
out the phonons-we give explicit expressions for these 
renormalized couplings below. 

The principal results of the present paper are Eqs. ifTjl 
and The latter is an analytic expression for the 

space-time spectral function that interpolates between 
the two limits in Eq. Q , and is shown to be very nearly 
exact in the exactly solvable case of forward scattering 
only. We later generalize it for el-ph couplings that in- 
clude backscattering interactions. Plots of A(k, E) com- 
puted from Eq. © are shown in the figures. 

There are several general features of A(k, E) that are 
worth noting. (1) As has been previously emphasized^ 
the spectral weight of the LL is concentrated in a roughly 
triangular region of the E-k plane (see Fig. 0), reflecting 
the fractionalized character of the elementary excitations, 
in contrast to the Fermi liquid case in which the spectral 
weight is concentrated along the line E = vp(k — kp), re- 
flecting the quasiparticle dispersion. (2) There is a renor- 
malization of the charge (holon) velocity produced by the 
el-ph coupling, such that v* < v c , which is analogous to 
the renormalization of Dp in the Fermi liquid. This pro- 
duces a "kink" in the spectrum, as shown in Fig. |S] (3) 
The E dependence of A(k, E) at fixed k is referred to 
in the ARPES literature as the energy distribution curve 
(EDC), while the k dependence at fixed E is referred 
to as the momentum distribution curve (MDC); the ex- 
tent to which there is a quasiparticle-like peak in the 
EDC of a LL is strongly dependent on the value of K c . 
For weakly interacting electrons, K c w 1, and the EDC 
exhibits a peak, although this peak contains a power- 
law tail indicating the absence of fermionic quasiparti- 
cles. For strong repulsive interactions that are sufficiently 
long range, K c <C 1, in which case the EDC is extremely 
broad, and can even fail to exhibit any well defined peak 
near the Fermi energy. In contrast, the structure of the 
MDC is much less variable^ and remains peaked even 
for small K c . The presence of extremely broad EDCs at 
the same time that the MDCs are narrow is a manifestly 
non-Fermi-liquid feature and a dramatic signature of the 
LL. This feature is illustrated in Fig. IT51 

For a model with purely forward scattering el-ph in- 
teractions, and for more general couplings as long as the 
cl-ph backscattering is not too strong, at low energies 
K c is increased by the presence of el-ph interactions, i.e. 
K* > K c . For the physically relevant case when the bare 
K c < 1, this means that features that are situated within 
luq of the Fermi energy are made sharper (more peaked) 



OJ, 
LU 

< 



no el-ph coupling (K c = 


0.15) A 


— — no el-ph coupling (K — 


0.3) '/ 




with el-ph coupling 

of strength ). = 0.75 


'/ 

1 1 




(K c = 0.15, K* = 0.3) 






i . i . i . i 


-,-s 

\ 

\ 

i . i 





-3 -2 

E/co„ 



FIG. 1: Comparison of the single-hole spectral functions at 
k — kF for a LL coupled to optical phonons (solid line), and 
a LL in absence of phonon coupling (dashed and dash-dotted 
lines). For the solid line, the el-ph coupling strength is A = 
0.75 (v c /v* — K*/K c — 2). For all curves, v c /v s = 4 and 
ujo/T — 10. See the plot for the values of K c . Normalizations 
are chosen for graphical clarity. 



by the coupling to phonons. Therefore, for K c <C 1, 
the result has the following similarity with the Fermi liq- 
uid case: features in A(k 7 E) appear broader for binding 
energies above ojq than for binding energies below. In 
the Fermi liquid case, this comes from a phonon-induced 
broadening at high energies, but in the LL it comes from 
a phonon-induced narrowing at low energies^ 

A plot of the EDC at k = kp is shown in Fig. ^ 
The solid line shows A(kp,E) for the case K c = 0.15, in 
the presence of forward scattering el-ph coupling, which 
produces, at low energies, a renormalized K* — 0.3. The 
dashed line shows A(kp,E) for the pure LL, but with 
K c = 0.3, while the dash-dotted line shows it for the 
pure LL with K c = 0.15. Clearly, for the LL coupled 
to phonons, the renormalized K c governs the properties 
at low binding energies, while the bare K c dictates the 
behavior at large binding energies. 

In Sec. [HJ we investigate the exactly solvable model 
(forward scattering only). The results are extended to 
general couplings using a renormalization group (RG) 
treatment in Sec. II I II In Sees. UVI andlVlwe speculate on 
the relevance of these results to real materials, especially 
the high-temperature superconductors. Appendix lAlcon- 
tains a derivation of the spectral function of the forward 
scattering only model. There we also present the exact 
result for the frequency- and momentum-dependent con- 
ductivity of this model. The optical conductivity is un- 
changed by the el-ph forward scattering, even though this 
interaction has dramatic effects on the spectral function. 
Appendix[B]contains technical details for the results pre- 
sented in Sec. IIIII 



3 



II. AN EXACTLY SOLVABLE MODEL 

The model is defined by the Hamiltonian density 

H = 7~Lll + 7~L P h + H e i- P h • (2) 
Here the purely electronic part of the Hamiltonian 



HhL = E v i 



K a Tl 2 a 



K a 



(3) 



is the famous spin-charge-separated Tomonaga-Luttinger 
liquid model of the interacting one-dimensional electron 
gas (1DEG) at incommensurate filling. It is expressed 
via bosonization, in terms of bosonic charge (a = c) and 
spin (a = s) fields <f> a , and their canonically conjugate 
momenta Tl a . Repulsive interactions usually renormal- 
ize the Luttinger parameter K c below its noninteracting 
value of 1 such that < K c < 1, and renormalize the 
velocities such that the v s < vp < v c . We will assume 
the system is spin-rotation invariant, which dictates that 
K s = 1. Expressions for the v Q 's and K a 's in terms of mi- 
croscopic short-range interaction parameters and reviews 
on the technique of bosonization can be found in many 
places in the Iiterature>i2ii2i22i2ii2£ Bosonization allows 
the fermionic fields to be expressed directly in terms of 
the bosonic fields: 



exp 



[rj((j) c + <T<fi s ) + C + , 



(4) 

where ^r),a is the right- or left-moving fermionic destruc- 
tion field (j] = ±1, respectively) for spin up or down 
(cr = ±1, respectively). Here a (x) = — J_ dx'U a (x') 
labels the dual bosonic field, and a is a short-distance 
cutoff corresponding to the lattice parameter. As is well 
known, the Hamiltonian in Eq. describes a line of 
fixed points for interacting electrons, and so captures the 
essential low-energy physics of a large class of physical 
systems. 

The purely vibrational part of the Hamiltonian 



W ph = [P 2 + ujIu 2 ] /2M 



(5) 



describes an Einstein oscillator. Here u and P are the 
phonon field and its canonical conjugate momentum, M 
is the ion mass, and again loq is the optical phonon fre- 
quency. Note that we work with units such that Boltz- 
mann's and Planck's constants are fc^ — Ti — 1. 

In the following section, we will consider a general form 
of the el-ph coupling, 7i i-ph, but for the purposes of the 
present section, we consider forward scattering interac- 
tions only, 



1-ph = OL 2 U p 



/ 2/7T U {d x (j)c 



(6) 



where a 2 is the el-ph coupling parameter, p is the long- 
wavelength component of the charge density, and the sec- 
ond equality makes use of the standard bosonization ex- 
pression for p. viet The forward scattering model is ex- 
actly solvable since in its bosonized form it is quadratic 



in the fields. In Appendix lAl we compute the renormal- 
ized couplings that define the {<?*} for this model. The 
exact result is 



where 



A= 2 



-KV c Mu! 2 



(7) 



(8) 



and the spin couplings are unrenormalized. Note that A 
depends on both the el-ph and el-el interactions. Also 
note that v c /v* = K*/K c > 1. At low energies, v c is 
reduced due to "phonon drag," while K c is increased 
due to an attractive interaction mediated by phonons. 
Since this model possesses an instability at A = 1, where 
the charge velocity goes to zero, A is restricted to the 
range < A < 1. The analog of this instability for the 
case of coupling to acoustic phonons has been studied 
previously^ 

In Appendix ^ we derive an exact expression for the 
spectral function. This quantity is most easily expressed 
in position space, but even then involves a momentum 
integral in the exponent that cannot be performed an- 
alytically. The integral can be evaluated numerically, 
and fortunately we are able to derive a simple analytic 
expression that accurately approximates it, as shown in 
Appendix lAl and Fig. We therefore use this analytic 
approximation in subsequent calculations. 

Specifically, our analytic approximation for the single- 
hole Green's function G v {x,t:X) = (^^(x, t)^ )ff (0, 0)) 
is 

G^tsAJ^G^tjO)^^!^, (9) 
g v {x,t; v c ,K c ,v c /ujo) 

where the exact Green's function in the absence of 
phonon coupling G v (x, t; 0) at temperature T = 1/0 is 24 

G v (x,t;0)= ^ Yl 9n{x,t;v a ,K a ,a) , (10) 



and we have defined the function 



g v (x,t;v,K, a) 

= TT (^sinh 
j=±i K 



ir(vt + jrjx — ia) 

v/3 



(11) 



Equation @ is the central result of this paper. The 
spectral function measured by ARPES is accurately given 
by its Fourier transform. 

Specifically, the single-hole spectral function for right- 
moving fermions is 25 



A(k,E) 



dx 



dte l(kx ' Et) dix^t-X) . (12) 
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FIG. 2: Dependence of the EDC at k — kF on A, 7 C , v c /v a , 
and ojo/T [(a), (b), (c), and (d) respectively]. In each panel, 
only one parameter is varied, and the rest are held fixed. 
Unless otherwise labeled, A = 0.75, j c = 0.6, v c /v s — 4, and 
uio/T — 10. For example, for the curves in (a), A has the 
value labeled in the plot legend, while y c = 0.6, v c /v s — 4, 
and uio/T — 10. In (a), the values of v c /v* = K* /K c are 1.69, 
2, 2.24 for A equal to 0.65, 0.75, 0.8, respectively. In (b), K c 
is 0.172, 0.150, 0.134 for 7 C equal to 0.5, 0.6, 0.7, respectively. 
Normalizations are chosen for graphical clarity (the apparent 
variation of the total spectral weight is an artifact). 



By combining this with Eq. |J3J, one obtains 

32 roo pec 

A(k, E) sa -r— / dx / dt X (x,T) 



TraujQ j_ oc J 
X cos [j3(kx-Et) -Q(x,i)] , (13) 

where we denned the dimensionless quantities [3 = 
ujoP/k, a = auJo/v c , 



k = v c (k - k F )/uj Q , E = E/uj , 
and the functions 

(ai/Vi) 2 



x(x,t)=n n 

i=ij=±i 

and 



sinh (i + jx/vi) + sin (di/vi) 



&(x,i) = 7ij arctan 

i=i j=±i 



tanh(i + jx/vi) 
t&n(ai/vi) 



with 



7ij = (1 - 2S i3 ) 



(Kj-3) 2 
8Ki 



(14) 



7y/2 



(15) 



(16) 



(17) 



K t = l, K 2 = K 3 =K C , 
«i = v s /v c , v 2 = v 3 = 1, 
ai = a,2 = a/ (3 1 a 3 = 04 



«4 = V*/V c , 



(18) 
(19) 
(20) 



Above, x and i are dimensionless dummy variables of 
integration, 5i 3 is the Kronecker delta, and a <^ 1 is a di- 
mensionless cutoff. The spectral function plots were ob- 
tained by performing the double integration in Eq. I|13l) 
numerically. Henceforth, we adopt the notation 



7c = (K c + K- 1 - 2)/8 = 721 



(21) 



which vanishes in the absence of el-el interactions. 

A(k, E) depends only on fc, E, loq/T, and on the three 
dimensionless parameters A, j c , and v c /v s m& In order to 
provide the reader with a qualitative understanding of 
how the spectral function depends on these parameters, 
in Fig. [3 we show how the EDC at k — kp changes when 
one is varied with the rest held fixed. Henceforth, we 
either chose representative parameters for the figures, or, 
for cases in which the theory is compared to experimental 
data, we fit the parameters. Not surprisingly, for cases 
in which we have fitted the parameters, they turn out to 
be somewhat material dependent. 

Figure[3Ja) exhibits the temperature dependence of the 
EDC at k = kp. The spectral weight of the low energy 
peak is reduced by increasing T, while at the same time 
the width of the peak increases proportional to T, due to 
the quantum critical nature of the LL. The value j c = 0.6 
(K c fa 0.15) used here indicates strong el-el repulsion. 

In Fig. 21 we show the dependence of the EDCs on fc, 
with and without forward scattering coupling to phonons, 
for uj /T = 10 and 7c = 0.6. Note the "double-peak" 
structure present near the phonon energy for moderate 
values of k. 

Figure [5] shows MDCs at various values of E for the 
same parameters as Fig. 0] A similar plot is shown in 
Fig. at a lower temperature (u>o/T = 40). Because 
of the lower temperature, here one can resolve three local 
maxima for moderate values of E. As shown in Fig. EJc), 
which is an enlargement of the E = — 1 curve, the cen- 
ter local maximum disperses at the renormalized charge 
velocity v*. Since this peak is by far the dominant one 
for \E\ <ti 1, the low-energy dispersion is characterized 
by what appears to be a single peak dispersing with ve- 
locity v*. At higher \E\, the v* peak disappears. If the 
temperature is increased sufficiently, the three local max- 
ima can no longer be resolved and appear as one peak, 
as seen in Fig. OJa). Also note that phonon coupling 
creates larger spectral weight at the spinon velocity v s , 
due to the higher effective K c (this feature disappears at 
high binding energies). 

In Fig. \7\ we present contour plots in the E-k plane, 
with and without phonon coupling, for the same param- 
eters as Fig. 0J Note the pronounced peak (red spot) at 
low binding energies in Fig. EJa), due to the increase in 
the effective K c . The reduction in charge velocity at low 
binding energies can also be seen here. 
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FIG. 3: Comparison of the temperature dependence of the 
EDCs at k = kp for (a) a LL coupled to phonons and (b) 
ARPES data of optimally doped Bi2212 (T c = 89 K) at the 
Fermi surface crossing along k == (0, 0) to (7r,7r) (the nodal 
direction) (Ref. Hij). In (a), ujq = 70 meV, 7 C = 0.6, and 
Vc/vs = 4. The solid line is for A = 0.75, and the dashed line 
is the result for the pure LL. 
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FIG. 5: MDCs (a) with and (b) without coupling to phonons 
of strength A = 0.72 for Wn/T = 10. All parameters are the 
same as Fig. H The curves are offset by E. 



Figure [8] shows the dispersion in the E-k plane, de- 
termined by fitting MDC curves to Lorentzian functions 
(the same method used in the ARPES literature). A 
change in the slope occurs at E w — 1. The ratio of the 
slope at \ E\ -C 1 to the slope at \E\ ^> 1 is approximately 
v*/v c = Vl — A. Note that the dispersion, when deter- 
mined by fitting to Lorentzians, is weakly T dependent. 



III. GENERAL ELECTRON-PHONON 
COUPLING 

In general, both forward and backward scattering (i.e. 
with momentum transfer near 2/cf) are possible. Thus, 
in general, we should consider both processes: 



K 



cl— ph 



«2 P + OL\ 



H.c 



(22) 



In this case, because $ is a nonlinear function of the 
bosonic fields, the problem is not exactly solvable. We 
therefore treat the backscattering term cti with a pertur- 
bative renormalization group scheme. It is important to 
note that if a\ is sufficiently strong, a gap opens up in the 
spin sector, and the system is a Luther-Emery liquid 2 ' 
(LEL) instead of a LL (see Appendix . We refer the 
reader to Ref. |2^ for a detailed study of the phase bound- 
ary separating the LL and LEL phases. Here we limit 
ourselves to the case in which the spectrum is gapless. 

It is convenient to define the dimensionless el-ph cou- 
plings 



FIG. 4: EDCs with (solid lines) and without (dashed lines) 
forward scattering coupling to phonons of strength A = 0.72 
(v c /v* = K*/K c = 1.9), with -y c = 0.6, v c /v s = 3, and 
ujo/T = 10. The vertical dashed line is drawn at the phonon 
energy. The notation is k — v c (k — fep)/wo and E = E/ujq, 
where E is measured with respect to the Fermi energy. The 
curves have vertical offsets proportional to k. 



At 



itvfMujq 



2 ' 



A 2 



itvfMluq 



2 ■ 



(23) 



In Appendix iBl we derive the following relations, which 
are generalizations of Eq. J7J) to the case of nonzero Ai: 



UcV /(l-A)(l-A + X c 2A), 



(24) 
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FIG. 6: MDCs (a) with and (b) without coupling to phonons 
for ujo/T = 40 and the same interaction strengths as Fig. 21 
(c) shows an enlargement of the E = — 1 curve from (a), and 
labels the velocities at which various features disperse. In (a) 
and (b) the curves are offset by E. 




FIG. 7: Contour plot for ujq/T — 10 of the single-hole spec- 
tral function (a) with and (b) without forward scattering to 
phonons of strength A = 0.72. All parameters are the same 
as in Fig. 31 



K* 



Kr 



1 - A 



1 - A + K?A 



(25) 



where we have introduced the effective backscattering el- 
ph parameter 



A = 



Vf 



2K r v, 



-A*, 



(26) 



which lies in the range < A < 1. Here A* is the renor- 
malized value of the bare el-ph backscattering parame- 
ter X\. For the case in which the el-ph interaction is 
unretarded (when the Fermi energy Ep < u>o), Ai re- 
mains unrenormalized (A* = Ai). However, for the phys- 




FIG. 8: Dispersion for A = 0.72, determined by fitting MDCs 
to Lorentzians, for uiq/T = 10 (solid line) and ujo/T — 40 
(dashed line). Interaction strengths are the same as in Fig. 
^] The dotted lines are fits to the high-|iJ| portion of the 
dispersions. 




FIG. 9: Dependence of the magnitude of the kink (v c /v*) 
and the degree to which the spectral function is peaked at 
low energies compared to high energies (K*/K c ) on the effec- 
tive backscattering el-ph interaction parameter A, for the case 
K c = 0.15. The dashed lines are the result in the absence of 
forward scattering el-ph interaction (A = 0), while the solid 
lines are for A = 0.75. 



ically interesting case Ep > ui , A* depends on the ratio 
Ep/ujQ and on the strength of the el-el interactions (see 
Appendix iBl for an explicit relation). As before, the for- 
ward scattering el-ph parameter is A = (2K c vf /v c ) A2; 
note that A2 contains no asterisk because it remains un- 
renormalized regardless of Ef/ojq^ 

Equations l|24() and (|25l) are actually completely gen- 
eral, nonperturbative expressions. However, our only 
way of relating the effective parameter A to the micro- 
scopic parameter Ai is through Eq. (|B6|) . which was ob- 
tained from one-loop RG, and is therefore valid only when 
the el-el interactions are weak and Ai <C 1. But Eqs. I|24|l 
and l|25(l remain valid even if the RG is carried out to an 
infinite number of loops. 

If fitting to particular ARPES data, the parameters 
K c , K*, and v*/v c are easily obtained. Then, A and 
A can be determined by inverting Eqs. (|24|l and (1231) . 
Therefore, for quasi-lD systems, ARPES is an effective 
probe of the relative amounts of forward and (renormal- 
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ized) backward scattering el-ph interactions. 

Note that for nonzero A, the relation K*/K c = v c /v* 
no longer holds; instead K*/K c = (1 — A)v c /v*. For 
K c < 1, the relation v* < v c holds regardless of A and A. 
However, K* > K c holds only if the ratio A/A is small 
enough. Specifically, if A/A > 1/(1+Jff), then if* < K c , 
which means that the low-energy spectral features are 
no longer made more peaked due to phonons. Figure [5| 
illustrates the result of increasing A at fixed A, for the 
case K c = 0.15. The magnitude of the kink, given by 
v c /v*, is increased by turning up A, while K* is reduced. 
For all the spectral function plots in this paper, we have 
set A = 0. 



IV. COMPARISONS WITH ARPES 
EXPERIMENTS IN THE CUPRATES 

In this section, we wish to illustrate the extent to which 
the observed ARPES spectra in the cuprate superconduc- 
tors resemble those of a LL coupled to optical phonons. 
At a gross level, the character of the resemblance be- 
tween the observed spectrum and a pure LL was estab- 
lished in Ref. 0. The strength of this analogy is further 
supported^ by direct comparison between the experi- 
mentally measured A(k, E) in a quasi-lD bronze and the 
cuprates. However, as more and better data have be- 
come available, it has become clear that there are fea- 
tures in the cuprate data-especially the widely reported 
"kinks" in the dispersion of the MDC peaka 1 ! 2 ! 30 ^ 1 ! 32 
that are qualitatively absent from the pure LL. The pure 
LL is also unable to reproduce the "nodal quasiparti- 
cles" seen in experiments while simultaneously producing 
broad features at high energy. Here, therefore, we pro- 
pose to make a comparison between the measured spec- 
tral functions, and the spectral function of a LL coupled 
to an optical phonon. 33 Since the LL describes a gapless 
state, we only compare our results to data taken in the 
gapless "nodal" direction, which is defined as k — (0, 0) 
to (7r,7r). In other directions, the ARPES spectrum of 
the cuprates develops a gap below a certain tempera- 
ture, whereas the spectrum remains gapless in the nodal 
direction even in the superconducting state. 

In Fig. 13 we have fitted the theoretical EDCs to 
ARPES data- in optimally doped Bi 2 Sr 2 CaCu208+,5 
(Bi2212) 34 at the Fermi surface crossing in the nodal di- 
rection, shown at various temperatures, both above and 
below T c . The experimental temperature dependence has 
been previously interpreted as evidence for quantum crit- 
ical behavior. Since the LL is a quantum critical state, 
the resemblance between theory and experiment in the 
present paper supports this interpretation. 

Figure ^| presents a fit to normal state EDCs for 
slightly underdoped 1 Bi2212 along the nodal direction, 
at various momenta. The theoretical curves contain a 
similar double-peak structure as the experiment. The 
line shapes of both the theory and experiment are char- 
acterized by extremely long high energy tails. 



(a) (b) 




-400 -40 

E (meV) E (meV) 

FIG. 10: Comparison of EDCs for (a) a LL coupled to 
phonons and (b) ARPES data of slightly underdoped Bi2212 
along the nodal direction (T c = 84 K) at various momenta 
(Ref. Q). For (a), A = 0.72, j c = 0.6, v c /v s = 3, uj = 70 meV, 
T = 80 K, and the dashed line is drawn at \E\ = uio- 




FIG. 11: Comparison of the dispersion of (a) the LL coupled 
to phonons with u>o = 70 meV at T = 30 K to (b) the disper- 
sion in Bi2212 (T c = 84 K) at T = 30 K in the nodal direction 
(Ref. U). For (a), all interaction strengths are the same as 
in Fig. 1101 For both plots, the dispersion was obtained by 
fitting MDC curves to Lorentzians. In (b), the authors define 
the rescaled momentum k' , by normalizing to 1 the value of 
kF ~ k at E = — 170 meV. The dashed lines are guides to the 
eye. 



Figure lllf a - ) shows the theoretical dispersion for the 
same parameters as Fig. llOf a). except with T = 
30 K, determined by fitting the theoretical MDCs to 
Lorentzians. Fig. Illf bl shows the experimental disper- 
sion for slightly underdoped 1 Bi2212 in the nodal direc- 
tion at T = 30 K, which the authors obtained with the 
same fitting procedure. The bizarre feature of the ex- 
perimental data in which the dispersion at high binding 
energies does not extrapolate to the origin is also seen in 
the theory^ 
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FIG. 12: Comparison of intensity plot of (a) the single-hole 
spectral function of a LL coupled to phonons with 7 C = 1.5, 
A = 0.93, ujq/T = 20, and v c /v s = 3 to (b) the ARPES 
spectrum of nonsuperconducting LSCO (3% doping) at T = 
30 K (Ref. 13(11) . In (a) and (b), the red line and points are the 
dispersion, obtained by fitting MDCs to Lorentzians. Panels 
(c) and (d) show enlargements of the dispersions; the dashed 
lines are fits to the high-binding-energy portion. 



For the theory plots, the manner by which the disper- 
sion is extracted (least-squares fitting to Lorentzians), 
yields a weighted average of the four velocities labeled in 
Fig. Efc). It is the contribution from the velocity —v c 
that is responsible for the plotted high-energy dispersion 
not extrapolating to the origin at E = and k = kp. 
This is because, when the dispersion is determined by fits 
to Lorentzians, the presence of spectral weight at k > kp 
[near the right side of the triangle in Fig. I12f a)] "pulls" 
the high- binding-energy dispersion to higher velocities. 

Figure I12f a) shows a contour plot of the theoreti- 
cal spectral function in the E-k plane and compares it 
to ARPES data2& for underdoped, nonsuperconducting 
La^Sr^CuCU (LSCO) with x = 0.03 [Fig. H2Tb)]. The 
value of A was chosen to give v c /v* = K*/K c = 3.8, 
which is the same as the ratio of the high to low binding 
energy velocities in the experimental plot. In Figs. 1121' c) 
and !12f d) we show the dispersions, obtained in the same 
manner as in Fig. ^] The EDCs and MDCs from these 
plots, at the Fermi momentum and Fermi energy, respec- 
tively, are shown explicitly in Fig. ^| For both the the- 
ory and experiment, the contrast between the sharpness 
of the MDCs and the breadth of the EDCs is dramatic. 

Note that the large values of j c and large values of 
A used for all the theoretical plots above indicate the 
presence of very strong el-el and el-ph interactions, with 
an effective el-ph interaction that is peaked in the forward 
scattering direction (since we used A = 0). 
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FIG. 13: EDC at k — kF for (a) LL coupled to phonons and 

(b) underdoped nonsuperconducting LSCO (Ref. 36) (nodal 
direction); and MDC at E — for (c) LL coupled to phonons 
and (d) underdoped LSCO (nodal direction). Panels (a) and 

(c) are slices from the plot in Fig. I12f a'): (b) and (d) are slices 
from the data in Fig. I12f b'). 



In Refs. |37] and |38|, an effort is made to fit the ex- 
perimental dispersion of the cuprates to the conventional 
theory of el-ph coupling in a Fermi liquid. Both papers 
report that getting a good fit to the experimental kink re- 
quires coupling between electrons and a broad spectrum 
of phonon modes that extends all the way up to two or 
three hundred meV. As the authors point out, this is 
unphysical, since from neutron scattering, the highest- 
energy phonon mode is less than 100 meV in all of the 
cuprates. This suggests that either phonons are not re- 
sponsible for the kink, or if they are, that the correct 
theory is very different from the conventional Fermi liq- 
uid one. In the present paper, the phonon density of 
states that couples to electrons is a single 5 function at 
u>o (Einstein phonon) . But even treating the phonon den- 
sity of states as adjustable, we suspect that a Fermi liq- 
uid treatment would not be able to fit the frequency and 
temperature dependence of the spectral function, nor the 
linear dependence of the MDC width on energyi 



V. CONCLUSIONS 

We have analyzed the effect of the el-ph coupling on the 
single-particle spectral function of the theoretically best 
understood (and exactly solvable) non-Fermi-liquid, the 
LL. Since, by definition, a non-Fcrmi-liquid is a state in 
which the elementary excitations are not simply dressed 
electrons or holes, A(k, E) should be the measurable 
quantity in which non-Fermi-liquid effects are most dra- 
matic. Thus, it is important to have a clear idea of which 
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features of this function best distinguish a Fermi liquid 
from a non-Fermi-liquid. It has been argued previously 
that the most dramatic signature of a LL is the appear- 
ance of extremely broad tails in the EDC, although peaks 
in the MDC remain relatively sharp. (This signature 
is particularly dramatic when the el-el interactions are 
strong; while the same distinction applies in principle for 
weak interactions, in practice the LL is harder to distin- 
guish from a Fermi liquid when the couplings are weak.) 

We have shown that some of the gross characteristic 
features of el-ph coupling in a Fermi liquid-an appar- 
ent kink in the dispersion relations and EDCs that are 
more peaked at \E\ < uj compared to \E\ > wo~are 
also present in a LL with strong el-el repulsion and for- 
ward scattering coupling to phonons. For a LL with more 
general el-ph couplings, the kink is still present, but the 
tendency of the EDC to be more peaked for \E\ < luq 
is eliminated by sufficiently strong el-ph backscattering. 
However, in either case, the Fermi liquid analogy can- 
not be taken too far: The basic discrepancy between the 
sharpness of the MDCs and the breadth of the EDCs, 
and the "triangular" confinement of spectral weight in 
the E-k plane, remain striking aspects of the LL which 
differentiate it from the Fermi liquid. 

We have also shown examples of measured spectral 
functions in the cuprates, and drawn attention to the 
similarities between them and our theoretical results. We 
believe that the similarities are dramatic. Exactly why 
the measured spectral functions look so much like the 
LL coupled to phonons is a deep question, which we will 
not explore further, here. However, at the very least, we 
feel that this comparison reinforces the conclusion, which 
has certainly been reached^&i on the basis of a variety 
of other experimental observations in the cuprates, that 
these materials are not well described in terms of the 
conventional electronic quasiparticles of simple metals. 
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a path integral over the bosonic fields n c , 4> c , Ii s , S , P, 
and it. Since the Lagrangian is quadratic in all of these 
fields, we can integrate out all fields except <j) c and <p s - 
This results in the Lagrangian C = C c [ip c ] + £ s [</? s ] with 



£ s VPs 



i 



2 2 \ 

v c q -A 

2 



2 2 2 



2K s v s 



H + vW] \<p,\ 



(Al) 
(A2) 



Here the bosonic Matsubara frequency is uj n — 2irn/fl, 
the field ip a = tp a (q,u> n ) is defined by <p a (x, r) = 
(27r/3L)- 1 E n /dg e-*"" T+i « x <p a (q,u n ), and L is the 
length of the system. The partition functional is then 
given by the path integral Z = J Vip c Vip s e - fi, <=[v<J-' s »IV.] 
with S a [ip a ] = (2nf3Ly 1 J2 n I d Q Calf a]- The zero- 
temperature dispersion relations for the hybridized 
charge-phonon bosonic collective modes are thus 



±V> 



v 2 q 2 ) 2 



2 v 2 q 2 



I 1/2 



(A3) 

Note that for A = 0, w_ = v c q is the long-wavelength 
density mode of the LL, and uj + — loq is the Einstein 
phonon mode. In Fig. I14f a1. we plot u>- and w+ for 
various values of A. 

It is a general feature of the bosonization approach that 
the fermionic correlation functions are most simply ex- 
pressed as a function of space and time. The Matsubara 
space-time Green's function was previously computed for 
a LL with forward scattering interactions with acoustic 
(instead of optical) phonons ^S4& The analytic structure 
of Eq. (|A1|I differs from the acoustic phonon model ac- 
cording to ujq — > cq, where c is the velocity of the acoustic 
phonon. The initial steps for computing the space-time 
Green's function involving functional integration over the 
fields and Matsubara frequency summations are the same 
for our model, except for this "parameter" change. We 
will therefore not reproduce these steps here, but rather 
refer the reader to Ref. 39. The final momentum inte- 
gral done in the exponent, however, is different for the 
two models. 

We write the single-particle imaginary-time Green's 
function for right moving fcrmions Q(x,t;X) — 

-^T T *i )ff (a;,r)*| )(T (0,0)^, where T T is the imaginary 

time ordering operator, as 



APPENDIX A: THE SPECTRAL FUNCTION 
AND CONDUCTIVITY OF THE EXACTLY 
SOLVABLE MODEL 

We now derive the spectral function of the model with 
purely forward scattering interactions. We work in the 
Matsubara representation, in which r = it is the imagi- 
nary time, and the partition functional is represented as 



2na 



exp[-/ c (5,f;A)-/ s (x,f)] (A4) 



(results for left-moving fermions are obtained by chang- 
ing x — ► —x). Here f c and f s are the charge and spin con- 
tributions to the exponent, which we will write in terms 
of the dimensionless variables x — xluq/v c and f = lu t. 
To simplify notation, for now we express the result in 
the limit T — > 0. For r > and infinitely large system 
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length, the exact result is 

/>oo 

f c (x,f;X) = / dqe- a « V{X[l-cos(p)e-^] 

(A5) 



i/=± 

iS w sin(p) e^" f }, 



f s (x,r) = F(x,t;v s /v c ,K s ,o), 
where a = aujQ/v c and we defined the functions 
(X c 2 + l)(c4-l) + A 



A± = 



B± = 



U)±-l 



2q{td 2 ± 



l + f± y/(l-q 2 ) 2 +4Aq 2 



(A6) 

(A7) 
(A8) 
(A9) 



and 



F(x,?;v,K,a) = ~ fx + 1 ] In 



x 2 + (a + vt) 2 



— — arctan 
2 



x 



a + vt 



(A10) 



For A = 0, the charge and spin parts are the same after 
an appropriate change of parameters: 

f c (x,f;0)=F{x,?;l,K c ,a), (All) 

yielding the known24 Green's function of a pure LL 



g(x,n o) 



2na 



Y[ H(x,T;v a ,K a ,a), (A12) 



where 

H(x, r; v, K, a) 



-i {K-iy/SK 



(a + vt) 2 + x 2 



a + vt — ix 



(A13) 



For arbitrary values of the parameters, although it is 
straightforward to evaluate Eq. (|A5|) numerically, it does 
not appear possible to perform the integration analyti- 
cally. We were, however, able to perform it analytically 
for the case of arbitrary A but a ^> 1: 



Q(x,t;X) 



2na 



H(x, t; v*,K*,a) H(x, r; v s , K s , a) 

for a > 1, (A14) 



with v* and K* given by Eq. 0. The limit a > 1 
is typically not satisfied in real materials, but since the 
spectral function is independent of a for |fc|, \E\ <C 1/a, 
we can use Eq. (|A14(I to deduce the exact behavior of 
the spectral function in the limit \k\, \E\ <C 1, regardless 
of the value of a: 

A(k, E; A) oc A LL (k, E; v*,K*,v s ,K s ) for \E\ « lj , 

(A15) 



where A^^k, E;v c , K c ,v s , K s ) = A(k,E;0) is the spec- 
tral function of the pure Luttinger liquid. 

For the physically interesting case a <C 1, we were 
able to derive an accurate analytic approximation for Eq. 
(|A5|I . which is the following: 

f c {x,f; A) w / c (x,f;0) 

+ F(x, f; v*/v c , K* c , 1) - F(x, f; 1, K c , 1), (A16) 
which gives for the Green's function 



G(x,t;X) ~G(x,t;0) 



H(x,t; v*,K*,v c /ujo) 
H{x,t;v c ,K c ,v c /uj q ) 

for a < 1, (A17) 



with v* and K* given again by Eq. Q. The T ^ 
single-hole real-time Green's function is written, using 
this approximation, in Eq. |(SJ. 

We demonstrate the accuracy of this analytic approx- 
imation by comparing Eq. (|A16|I with the exact result, 
obtained by performing the integration in Eq. (|A5I) nu- 
merically. From Figs. I15f c) and I15f d^ we see that for 
A = 0.75 and 7 C = 0.6, the approximation gives within 
9% of the exact result for the real part of f c (x, f; A). For 
A = 0.25 [Figs. I15f a) and I15f b1]. the error is less than 
3%. The agreement with the imaginary part is similarly 
excellent. 

We have also computed the frequency- and 
momentum-dependent conductivity for the LL with 
forward scattering off phonons. The exact result for the 
real part of the conductivity at T = is 



cri(q,w;A) = e 2 K c v c }^ [WS(uj - vu}+) 

L/=±l 

+ (l-W)5{u-vwJj\. 



(A18) 



Here e is the charge of the electron, 8 is the delta function, 
and the spectral weight of the w+ mode is 



W 



(A19) 



which depends only on v c q/uJo and A. We plot W in 
Fig. I14f b). Note that the conductivity at q = (optical 
conductivity) remains unchanged by forward scattering 
off phonons: 



£Ti(Q,w;A)=(ri(0,w;0) = 2e 2 K c v c 8{u), 



(A20) 



yet this same interaction produces dramatic changes to 
the spectral function. 



APPENDIX B: RG TREATMENT OF GENERAL 
ELECTRON-PHONON COUPLING 

In the case in which the el-el couplings and the 
backscattering el-ph coupling Ai [Eq. I)23|l] are weak, 
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FIG. 14: Panel (a) shows the dispersion of the bosonic modes 
w+ (upper curves) and ui- (lower curves), (b) shows the spec- 
tral weight W of the uj+ mode in the conductivity. The nota- 
tion is lu± = u)±/gjo and q — v c q/uo. 



X = 0.25, x= 1 



• exact 
■ approximate 



X = 0.25, x= 10 




FIG. 15: Comparison of the exact Green's function charge 
exponent f c (x,f; A), obtained by evaluating Eq. <A5t nu- 
merically (solid lines), with the approximate analytic expres- 
sion in Eq. (IA16I (dashed lines). Panels (a) and (b) are for 
A = 0.25; (c) and (d) are for A = 0.75. For all panels, 7c = 0.6 
and a = 0.05. Panels (a) and (c) show the f dependence for 
5 = 1; (b) and (d) are for x = 10. The main plots show the 
real part of f c and the insets show the imaginary part. For 
the imaginary part, the dashed and solid lines overlap. 



although of arbitrary relative strength, the phonon- 
induced renormalization of the LL parameters can be 
computed using a two-step perturbative (one-loop) RG 
scheme Since, for Ep 3> ^o, the coupling Ai is 
strongly renormalized, and the RG flows of Ai are modi- 
fied by the direct el-el interactions, the resulting expres- 
sions for the renormalizations of v c and K c are more 
complicated than in Eq. Q. However, they have re- 
cently been analyzed in detail by one of uSf2& and will be 
summarized below. 

We employ standard notation for the important short- 
range el-el interaction parameters of the incommensu- 
rate 1DEG: gi (backscattering), g 2 (forward scattering 
on both left- and right-moving branches), and g^ (for- 
ward scattering on only one branch). The Hamiltonian 
density for the el-el interaction portion is then 



H e l-el = gi ^1^-1,^1,^-1,^ 
<r,a'=±l 

+ 92 J2 *l,«r*il lff '*-l,«r'*l,«r 
<r,a'=±l 

+ 34 E 

Vt (T= it 



vtrt \[/t m m 



(Bl) 



For the extended Hubbard model in the continuum 
(weak-coupling) limit, these parameters are g\ = U— 2V', 
g 2 = U + 2V', and g 4 = U/2 + 2V', where U is the on- 
site interaction, V' = — V cos(2kp), and V is the nearest 
neighbor interaction (near halffilling, V' sw V). 

For an incommensurate 1DEG in the presence of a gen- 
eral el-ph coupling, the effective charge parameters at low 
energy are 



K* 



v F ^{l+grf-{g^/2f 



y i + 5i ot -5^ ot /2' 

where we defined 

- gi — 2g 2 — (AJ - 2A 2 



9c 

tot 

9i 



94 - M 



(B2) 
(B3) 



(B4) 
(B5) 



and g~i — gi/irvp. Here the renormalized el-ph backscat- 
tering coupling is denoted by A*. For the case Ep ^> ljo 
and gi, g 2 , Ai <C 1, it is given bj^ 



K = 



Xih(l ) 



where 



h{x) 



1 — Ai J ° dx h(x) 

exp[-(g! - 2g 2 )x/2] 
(1 + g lX ) 3 / 2 



(B6) 



(B7) 



Iq = \n{Ep /ujq), and Ai and A2 are the bare el-ph cou- 
plings defined in Eq. I|23|) . 
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In the weak-coupling limit and for Ep luq, the total 
effective backscattcring interaction after integrating out 
degrees of freedom from Ep to loq is given by 

g r = -^—-X\. (B8) 

In order for the Luttinger liquid description to be valid, 
it is required that g\ ot is repulsive (g\ ot > 0). In this 
case, the effective spin parameters are K* = 1 and v* — 
vp{l — gi). However, if A* is large enough to cause g[ ot < 
0, the RG flows carry K* and v* to zero at low energies, 
signaling the emergence of a spin gap. In this case, the 
system is described as a Luther-Emery liquic— Z phase. 



We refer the reader to Ref. [2^ for a detailed study of the 
phase diagram of the 1DEG coupled to phonons. Since 
the spectral function computed in the present paper is 
not applicable to a Luther-Emery liquid, here we restrict 
our attention to the case in which there is no spin gap. 

By combining the above expressions for v* and K* with 
the expressions for v c and K c [given by Eqs. (|B2|I and 
(|B3|) with A* = A2 = 0], one can rewrite the results in a 
form that shows more clearly how the relation v*/v c — 
K c /K* = Vl — A is modified by the presence of el-ph 
backscattering-see Eqs. (|24() and l|25[) . There, the effects 
of el-ph backscattering are contained in the parameter A, 
defined in Eq. J2BJ. 
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